Late-Time Tails of Wave Propagation in Higher Dimensional Spacetimes 
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We study the late-time tails appearing in the propagation of massless fields (scalar, electromag- 
netic and gravitational) in the vicinities of a D-dimensional Schwarzschild black hole. We find that 
at late times the fields always exhibit a power-law falloff, but the power-law is highly sensitive to 
the dimensionality of the spacetime. Accordingly, for odd D > 3 we find that the field behaves as 
^-{2i+D-2) jg^^g times, where I is the angular index determining the angular dependence of the 
field. This behavior is entirely due to D being odd, it does not depend on the presence of a black 
hole in the spacetime. Indeed this tail is already present in the flat space Green's function. On the 
other hand, for even D > 4 the fleld decays as ^-(^i+SD-s)^ ^^^^ ^jjjg time there is no contribution 
from the fiat background. This power-law is entirely due to the presence of the black hole. The 
D = 4 case is special and exhibits, as is well known, the ^~(2i+3) \yg.\iQYioi. In the extra dimensional 
scenario for our Universe, our results are strictly correct if the extra dimensions are infinite, but 
also give a good description of the late time behaviour of any field if the large extra dimensions are 
large enough. 

PACS numbers: 04.70.-s, 04.30.Nk, 04.70.Bw, ll.25.-w 



I. INTRODUCTION 



It is an everyday life experience that light rays and 
waves in general travel along a null cone. For example, if 
one lights a candle or a lighter for five seconds and then 
turns it off, any observer (at rest relative to the object) 
will see the light for exactly five seconds and then sud- 
denly fade out completely. Mathematically this is due 
to the well known fact that the flat space 4-dimensional 
Green's function has a delta function character and there- 
fore has support only on the light cone. There are how- 
ever situations where this only-on-the light cone prop- 
agation is lost. For instance, in a curved spacetime a 
propagating wave leaves a "tail" behind, as shown by 
DeWitt and Brehme's seminal work [1]. This means that 
a pulse of gravitational waves (or any massless field for 
that matter) travels not only along the light cone but 
also spreads out behind it, and slowly dies off in tails. 
Put it another way, even after the candle is turned off in 
a curved spacetime, one will still see its shinning light, 
slowly fading away, but never completely. This is due to 
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backscattering off the potential [2] at very large spatial 
distances. 

The existence of late-time tails in black hole spacetimes 
is by now well established, both analytically and numeri- 
cally, in linearized perturbations and even in a non-linear 
evolution, for massless or massive fields [2-9]. This is a 
problem of more than academic interest: one knows that 
a black hole radiates away everything that it can, by the 
so called no hair theorem (see [10] for a nice review), but 
how does this hair loss proceed dynamically? A more or 
less complete picture is now available. The study of a 
fairly general class of initial data evolution shows that 
the signal can roughly be divided in three parts: (i) the 
first part is the prompt response, at very early times, and 
the form depends strongly on the initial conditions. This 
is the most intuitive phase, being the obvious counter- 
part of the light cone propagation, (ii) at intermediate 
times the signal is dominated by an exponentially decay- 
ing ringing phase, and its form depends entirely on the 
black hole characteristics, through its associated quasi- 
normal modes [11, 12]. (iii) a late-time tail, usually a 
power law falloff of the field. This power law seems to 
be highly independent of the initial data, and seems to 
persist even if there is no black hole horizon. In fact it 
depends only on the asymptotic far region. Mathemati- 
cally each of these stages has been associated as arising 
from different contributions to the Green's function. The 
late-time tail is due to a branch cut [5]. The study of 
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linearized (wc note that non-linear numerical evolution 
also displays these tails, but here we shall work at the 
linearized level) perturbations in the black hole exterior 
can usually be reduced to the simple equation 

[d1-dl + V{x)]^ = Q, (1) 

where the potential V{x) depends on what kind of field 
one is considering and also, of course, on the spacetime. 
A detailed study of the branch cut contribution by Ching, 
Leung, Suen and Young [6, 7] has provided analytical re- 
sults for some specific but quite broad class of potentials. 
These analytical results concerning the late time tails 
were confirmed numerically. 

It is not generally appreciated that there is another 
case in which wave propagation develops tails: wave 
propagation in odd dimensional flat spacetimes. In fact, 
the Green's function in a i?-dimensional spacetime [IS- 
IS] have a completely different structure depending on 
whether D is even or odd. For even D it still has sup- 
port only on the light cone, but for odd D the support of 
the Green's function extends to the interior of the light 
cone, and leads to the appearance of tails. It is hard 
to find good literature on this subject, but a complete 
and pedagogical discussion of tails in flat D-dimensional 
backgrounds can be fomid in [15]. 

A study of wave physics in higher dimensions is now, 
more than ever, needed. It seems impossible to for- 
mulate in four dimensions a consistent theory which 
unifies gravity with the other forces in nature. Thus, 
most efforts in this direction have considered a higher 
dimensional arena for our universe, one example being 
string theories which have recently made some remark- 
able achievements. Moreover, recent investigations [16] 
propose the existence of extra dimensions in our Uni- 
verse in order to solve the hierarchy problem, i.e., the 
huge difFcrcncc between the electrowcak and the Planck 
scale, TOEw/A'/pi ^ 10"^''. The fields of standard model 
would inhabit a 4-dimensional sub-manifold, the brane, 
whereas the gravitational degrees of freedom would prop- 
agate throughout all dimensions. 

A first step towards the understanding of gravitational 
wave physics in higher dimensions was given by Cardoso, 
Dias and Lcmos [14], by studying wave generation and 
propagation in generic £)-dimensional flat spacetimes. 
Here we shall take a step further, by studying wave tails 
in higher dimensional black hole spacetimes. We will re- 
strict the analysis to higher dimensional Schwarzschild 
black holes. As expected, if one now considers tails in 
higher dimensional black hole spacetimes two aspects 
should emerge: in odd dimensional spacetimes one ex- 
pects the black hole contribution to the tail to be smaller 
than that of the background itself. Therefore for odd di- 
mensions the tail should basically be due to the flat space 
Green's fimction. However, for even D-dimcnsional black 
hole spacetimes there is no background contribution, and 
one expects to see only the black hole contribution to the 
tail. A recent study by Barvinsky and Solodukhin [17] 
has showed that such tails may not be impossible to de- 



tect. Unfortunately, the weakness of gravitational waves 
impinging on Earth, make this an unlikely event. We 
note however that they worked with small length, com- 
pact extra-dimensions, whereas we shall consider large 
extra dimensions. Our results will be strictly correct if 
the extra dimensions are infinite, but also allow us to de- 
termine the correct answer if the large extra dimensions 
are large enough that the timescale for wave reflection at 
the boundaries is larger than the timescales at which the 
tail begins to dominate. 

The evolution problem in a i3-dimensional 
Schwarzschild background can be cast in the form 
(1), and we will show also that the potential can be 
worked out in such a way as to belong to the class of 
potentials studied in [6, 7]. Therefore, their analytical 
results carry over to the D-dimensional Schwarzschild 
black holes as well. We will verify this by a direct 
numerical evolution. The main results are: the late-time 
behavior is dominated by a tail, and this is a power-law 
falloff. For odd dimensions the power-law is determined 
not by the presence of the black hole, but by the fact 
that the spacetime is odd dimensional. In this case the 
field decays as ^-(2;+d-2)^ where I is the angular 

index determining the angular dependence of the field. 
This is one of the most interesting results obtained 
here. One can show directly from the flat space Green's 
function that such a power-law is indeed expected in flat 
odd dimensional spacetimes, and it is amusing to note 
that the same conclusion can be reached directly from 
the analysis of [6, 7]. For even dimensional spacetimes 
we find also a power-law decay at late times, but with 
a much more rapid decay, \1/ ~ ^-(2i-i-3D-8) _ -poY even 
D, this power-law tail is entirely due to the black hole, 
as opposed to the situation in odd D. These results are 
strictly valid for D > 4. Four dimensional Schwarzschild 
geometry is special, having the well known power-law 
tail * ~t-(2i+3). 



II. A BRIEF SUMMARY OF PREVIOUS 
ANALYTICAL RESULTS FOR A SPECIFIC 
CLASS OF POTENTIALS 

In a very complete analysis, Ching, Leung, Suen and 

Young [6, 7] have studied the late-time tails appearing 
when one deals with evolution equations of the form (1), 
and the potential V is of the form 

T./ N I'il' + l) , Ci\0gX + C2 , , 

^(^) ^ + ^ ,a;^oo. (2) 

By a careful study of the branch cut contribution to the 
associated Green's function they concluded that in gen- 
eral the late-time behavior is dictated by a power-law or 
by a power-law times a logarithm, and the exponents of 
the power-law depend on the leading term at very large 
spatial distances. The case of interest for us here, as we 
shall verify in the following section, is when ci = 0. Their 
conclusions, which we will therefore restrict to the ci = 
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case, are (see Table 1 in [6] or [7] 



(i) if V is an integer the term \^ does not contribute where S1_d-2 



to the late-time tail. We note this term represents just 

the pure centrifugal barrier, characteristic of flat space, so 
one can expect that indeed it does not contribute, at least 
in four-dimensional spacetime. We also note that since 
even dimensional spacetimcs have on-light cone propaga- 
tion, one may expect to reduce the evolution equation to 
a form containing the term '^^^t"'"'' with p an integer. We 
shall find this is indeed the case. Therefore, for integer u, 
it is the term that contributes to the late-time tail. In 
this case, the authors of [6, 7] find that the tail is given 
by a power-law, 

^ ^t~^ , fj.>2u + a , aodd integer <2i' + 3. (3) 

For this case (a an odd integer smaller than + 3) the 
exponent fj, was not determined analytically. However, 
they argue both analytically and numerically, that fx = 
2v + 2a — 2. For all other real a, the tail is 



^ ~ i-Ca-'-K") ^ all other real a. 



(4) 



(ii) if V is not an integer, then the main contribution to 
the late-time t; 
case the tail is 



the late-time tail comes from the '^'^'^t^'' term. In this 



non — integer v. 



(5) 



We will now see that for a £)-dimensional Schwarzschild 
geometry the potential entering the evolution equations 
is asymptotically of the form (2) and therefore the results 
(3)-(5) can be used. 



III. THE EVOLUTION EQUATIONS AND 
LATE-TIME TAILS IN THE D-DIMENSIONAL 
SCHWARZSCHILD GEOMETRY 

Here, we shall consider the equations describing the 
evolution of scalar, electromagnetic and gravitational 
weak fields outside the D-dimensional Schwarzschild ge- 
ometry. We shall then, based on the results presented in 
the previous section, derive the late-time tails form of the 
waves. We will find they are always a power-law falloff. 



A. The evolution equations and the reduction of 
the potential to the standard form 

The metric of the D-dimensional Schwarzschild black 
hole in {t,r,9i,92, ..,9d-2) coordinates is [18] 



with 



/ = 1-^ 

■' rD-3 



(6) 
(7) 



The mass of the black hole is given by ?6^s 

is the area of a unit {D — 2) 



^(0-l)/2 



■ r[(D-i)/2] 

sphere, and dfljj_2 is the line element on the unit sphere 
5^-2 ■\^g only consider the linearized approxima- 
tion, which means that we are considering wave fields 
outside this geometry that are so weak they do not alter 
this background. Technically this means that all covari- 
ant derivatives are taken with respect to the background 
metric (6). The evolution equation for a masslcss scalar 
field follows directly from the (relativistic) Klein-Gordon 
equation. After a separation of the angular variables with 
Gegenbauer functions (see [19] for details) we get that the 
scalar field follows (1) with a potential 



Vs{r.) = f{r) 



(D-2)(£)-4)/(r) , {D^2)f'{r) 



4^2 



-I- 



2r 



(8)" 

where r is a function of the tortoise coordinate r^, accord- 
ing to ^ = /(r). The constant a = Z(/ + D — 3) is the 
eigenvalue of the Laplacian on the hypersphere S^~^, 
and /'(r) = I can take any nonnegative integer 

value. Of course the evolution equation is (1) where the 
variable x is in this case the tortoise coordinate r, . This 
is the standard form in which the potential is presented. 
However, one can collect the different powers of r and get 



V.{r.) = f{r) 



v{v + 1) \ {D- 2fM 



+ 



where 



(9) 



(10) 



Asymptotically for large r* one can show that 



K(n) 



v{v + 1) I [D- 2)Ml 



+ 



{3-1- D). 



£>-4 

(11) 

This is strictly valid for D > 4. In the D = 4 case 

there is a logarithm term [2]. Notice that the coefficient 
appearing in the centrifugal barrier term '^'■^j"^-* is, as 
promised, an integer for even D, and a half-integer for 
odd D. The gravitational evolution equations have re- 
cently been derived by Kodama and Ishibashi [20]. There 
are three kinds of gravitational perturbations, according 
to Kodama and Ishibashi's terminology: the scalar grav- 
itational, the vector gravitational and the tensor gravi- 
tational perturbations. The first two already have their 
counterparts in _D = 4, which were first derived by Regge 
and Wheeler [21] and by Zerilli [22]. The tensor type is 
a new kind appearing in higher dimensions. However, it 
obeys exactly the same equation as massless scalar fields, 
so the previous result (9)-(ll) holds. It can be shown in 
fact that the scalar and vector type also obey the; same 
evolution equation with a potential that also has the form 
(11) with a slightly different coefficient for the l/rf~^ 
term. For example, for the vector type the potential is 

'v{v + l) 1 3{D-2)^M' 



VM = fir) 



„D-l 



(12) 
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where v is defined in (10). Therefore asymptoticaUy for 
large , 



Fgv(r*) 



v{v + \) 1 {D-2)M 



„D-1 



[8 + £'2 + D(Z-6)+;(Z-3)] 



(13) 



which is of the same form as the scalar field poten- 
tial. The scalar gravitational potential has a more com- 
plex form, but one can show that asymptotically it has 
again the form (11) or (13) with a different coefficient in 
the 1/rf ~^ term. Since the explicit form of this coeffi- 
cient is not important here, we shall not give it explic- 
itly. Electromagnetic perturbations in higher dimensions 
were considered in [23]. Again, asymptotically for large 
r* they can be reduced to the form (11) (where again 
v = I — 2 + Y and the 1/rf ~^ coefficient is different) , so 
we shall not dwell on them explicitly. 



B. Late-time tails 

Now that we have shown that the potentials appear- 
ing in evolution of masslcss fields in the D-dimensional 
Schwarzschild geometry belong to the class of potentials 
studied in [6, 7] we can easily find the form of the late- 
time tail. For odd dimensional spacetimes, v is not an 
integer, therefore the centrifugal barrier gives the most 
important contribution to the late-time tail. According 
to (5) we have that the late-time tail is described by the 
power-law falloff 



^^t-(2'+0-2),oddD. 



(14) 



According to the discussion in the introduction, this tail 
is independent of the presence of the black hole, and 
should therefore already appear in the flat space Gkkui's 
function. Indeed it does [14, 15]. The flat, odd dimen- 
sional Green's function has a tail term [14, 15] propor- 
whcrc 9 is the Heavisidc step 



e(t-r) 



tional to ^^2_r^)n/2-l ■ 
function. It is therefore immediate to conclude that, 
for spherical perturbations, for example, the tail at very 
large times should hct^'-^^^^ which is in agreement with 
(14) for / = (which are the spherical perturbations). It 
is amusing to note that the analysis of [6, 7] gives the 
correct behavior at once, simply by looking at the cen- 
trifugal barrier! We have checked numerically the result 
(14) for D = 5, and the results are shown in Figs. 1 and 
2. 

The numerical procedure followed the one outlined in 

[3], with constant data on v = Vq, where v is the ad- 
vanced time coordinate, v = t + r^,. One final remark 
is in order here. When numerically evolving the fields 
using the scheme in [3], we have found that for _D > 5 
the tail looked always like t~^^'^~^^\ This is a fake be- 
havior, and as pointed out in [7] (see in particular their 
Appendix A) it is entirely due to the ghost potential, 
appearing for potentials vanishing faster than for a 
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FIG. 1: Generic time dependence of a scalar field ^ in a 
five dimensional Schwarzschild geometry, at a fixed spatial 
position. We took as initial conditions a Gaussian wave packet 
with (7 = 3 and Vc = 10. We have performed other numerical 
extractions for different initial values. The results for the 
late-time behavior are independent of the initial data, as far 
as we can tell. For I = the late-time behavior is a power- 
law with * t~^-^, for / = 1, ~ t~^-^ at late times and for 
Z = 2, ^I' ~ t"'^'^. The predicted powers are —3, —5 and —7, 
respectively. 
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FIG. 2: Generic time dependence of a gravitational Gaussian 
wavepacket in a five dimensional Schwarzschild geometry, 
at a fixed spatial position. The results for the late-time be- 
havior are independent of the initial data. For I = 2 the 
late- time behavior is a power-law with ^' ~ for Z = 3 

the falloff is given by ^' ~ t~^ '^ at late times, and for Z = 4 
it is ^' ~ t~^^'^. The predicted powers are —7, —9 and —11, 
respectively. 
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FIG. 3: Evolution of a field subjected to a fivo-dimensional 
model potential. The potential is given by (9) but here we 
take as an integer, v = 1. We obtain at late times ^' ~ t~^'^, 
whereas the analytical prediction is t~^. This shows nicely 
that for integer u the centrifugal barrier contribution vanishes, 
and it is the next term in the asymptotic expansion of the 
potential that gives the most important contribution. We 
have checked that this power law is indeed the correct one, 
and not a numerical artifact of the ghost potential. 

second order scheme. Technically the presence of ghost 
potentials can be detected by changing the grid size [7]. 
If the results with different grid sizes are different, then 
the ghost potential is present. Our numerical results for 
D — b, presented in Figs. 1-3 are free from any ghosts. 

The numerical results are in excellent agreement with 
the analytical predictions (14), and seems moreover to 
be quite independent of the initial data. This also means 
that for odd dimensional spacetimes the late time behav- 
ior is dictated not by the black hole, but by the fact that 
spacetime has an odd number of dimensions. To further 
check that it is in fact the centrifugal barrier term that 
is controlling the tail, we have performed numerical evo- 
lutions with a five dimensional model potential. To be 
concrete, we have evolved a field subjected to the poten- 
tial (9) (with = 5) but we have considered an integer 
value for v, namely v = 1. The result is shown in Fig. 
3. Of course the true potential has a semi-integer value 
for V, but this way one can verify the dependence of the 



tail on the centrifugal term. Indeed ii v = \ then by 
(4), with 11 = 1 and a = 4, the late time tail should be 
^ ~ . The agreement with the numerical evolution is 
great. It is therefore the centrifugal barrier that controls 
the tail in odd dimensional spacetimes. For even dimen- 
sional spacetimes, v — I — 2 + !^ is an integer. Moreover, 
a = D — 1 < 21/ + 3 = 21 + D — 1. Therefore we are in 
situation (3). So the late time tail of wave propagation 
in an even D-dimensional Schwarzschild spacetime is a 
power-law, 

*~f-(2'+3i^-8),even£). (15) 

IV. CONCLUSIONS 

We have determined the late time behavior of mass- 
less fields (scalar, electromagnetic and gravitational) out- 
side a D-dimensional Schwarzscliild black hole. For 
odd D, the field at late times has a power-law falloff, 

~ ^-(2(-i-D-2)^ ^j^jg ^^jj jg independent of the pres- 
ence of the black hole. It depends solely on the flat space- 
time background, through the properties of the flat space 
odd dimensional Green's function. For even D, the late 
time behavior is again a power-law but this time it is 
due to the presence of the black hole, and is given by 

^ f-i2i+3D-8)^ Ig^^g ^.jjj^gg fQj. £) > 4. -yVe have 

focused on large extra dimensions only. Recent investi- 
gations [17], focusing on brane world models, and there- 
fore compact extra dimensions suggest that the tail is 
more slowly damped if the extra dimensions arc compact. 
In fact, for £) = 5 they obtain a power-law \1/ ~ t"^/^ 
whereas we have ^ ~ for spherically symmetric per- 
turbations. This may be due to the reflection of the field 
at the boundaries of the extra dimension. 
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